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Abstract 


The analytical solution for the Stokes vector of a light field inside and 
outside of a disperse asymmetric isotropic media is obtained. It is assumed 
that scatterers are contained in a plane-parallel vertically and horizontally 
homogeneous slab. The diffused light field is studied in the framework of 
the single-scattering approximation. Equations obtained can be applied for 
the development of radiative transfer codes for isotropic asymmetric media, 


based on the adding—doubling procedure. 
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1. Introduction 


Polarized light transport in random media is usually studied 
in the framework of the vector radiative transfer theory. 
This theory is based on solutions of the system of integro- 
differential equations, which can be written in the following 
form [1, 2]: 


dS(Q) 
dx 


zega Ge f dÈ êsa BSG’) O) 


An 


where the vector Q describes the direction of li ght propagation, 
S is the Stokes vector, Ĝext and Gycq are extinction and 
differential scattering matrices and x is the distance. It is 
assumed here that both medium and light source are stationary 
and there are no internal sources of radiation inside the 
medium. Also, the possible change of frequency in the 
scattering process is neglected. Equation (1) follows from 
the Maxwell theory, if one allows for several simplifying 
assumptions [3]. 

Equation (1) can be applied to disperse media of arbitrary 
shapes. Experimentalists, however, deal mostly with plane- 
parallel slabs. Thus, only a transport problem in the 
plane-parallel geometry is studied here. We assume that a 
random medium (e.g. a suspension of macroscopic particles) 
is contained in a plane-parallel slab, which extends to infinity 
in the horizontal plane. However, a finite thickness along the 
vertical coordinate axis O Z is assumed (see figure 1). Then itis 
convenient to measure distances in the direction perpendicular 
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to the border of the medium and x = z/cos V, where V is the 
zenith observation angle, counted from the downward vertical, 
and z is the Cartesian coordinate along axis OZ. The vector 
Qin equation (1) in this case is determined by the zenith angle 
v and azimuth angle @. 

We will assume that the medium is uniformly illuminated 
from the top by a infinitely wide light beam of an arbitrary 
polarization. The angle of incidence, measured from the 
vertical, is equal to Vo. Then equation (1) is written in the 
following form [2]: 


dS, $) 
Os 0 ————— 
dz 


f a f do’ sin ¥ Gsca(d, O, 6, PDS, p. (2) 
0 0 


= Cox S (Ù, H) 


The task of this paper is to find the solution of equation (2), 
assuming that secondary and higher-order scattering events 
can be neglected. Also, we make the assumption of the 
homogeneity of the random medium under consideration. This 
means that extinction and scattering matrices are constant 
throughout the scattering layer studied. 

A similar problem was considered by Ishimaru and 
Yeh [4]. However, they provided the solution only for a light 
field escaping from the medium. We are interested in the light 
field inside the slab as well. 

We also study the special case of turbid media with 
optically active spherical particles. 
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Figure 1. The geometry of the problem. 


2. Coherent beam 


The Stokes vector Š (Q) in equation (2) can be presented as the 
sum of two parts: 


S(Q) = TQ) + 1,.(Q)5(Q — Qo) (3) 


where the Stokes vector I (Q) describes the incoherent or 
diffused light beam and AOO — Qo) is the Stokes vector 
of the coherent or direct beam propagated in the direction 
Q0(%, po). Here 8(Q — Qo) is the delta function. 

Let us substitute equation (3) into (2). Then we have two 
separate equations: one for the diffused light beam Stokes 
vector T (V, ġ) and one for the direct light beam Stokes vector 
I.0, $): 


OS p Oo) = Bex I (8, o) 


dz 
š f ag | dd’ sin P'ôsa(, D, o, TOG) 
0 0 
+ sca (o, Po, D, O)T: (O0, Go) (4) 
and 7 
d7. (vo, n > 
awaa E tr N (5) 


dz 
Equation (5) can be easily solved [4]. The solution obtained 
should be substituted into equation (4). In particular, we obtain 
from equation (5): 7 7 
L=ÎJ, (6) 
where 
T = exp(—Gexil), (7) 


l = z/ cos Yo and J is the Stokes vector of the incident light. 
Equation (4) can be written in the following form to 
account for equation (6): 


os P $) 


a —Gexr1 O, p) 
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2n us 
+ / ay’ f d’ sin Y 6sca(d, O, 6, POI, 6’) 
0 0 


+ Bsca(Bo, Po, V, PIT I (Vo, Qo). (8) 


As far as boundary conditions are concerned we assume that 
there is no diffused light entering a scattering layer from 
outside. 

Equation (8) is a system of four coupled integro- 
differential equations. They can be used for studies of 
polarized diffused light transfer both in symmetric and 
asymmetric media. For symmetric media, the centre of any 
large spherical volume is a centre of symmetry and any plane 
through this centre is a plane of symmetry [5]. This is not 
the case for asymmetrical media, which usually have some 
rotatory power. We assume in this study that an asymmetrical 
medium is macroscopically isotropic, however. The isotropy 
could in principle be due to an isotropic distribution of small 
anisotropic elements. 

Clearly, equation (8) can be solved using only numerical 
techniques. A good starting point for the numerical procedure 
would be the solution of equation (8) in the single scattering 
approximation [4], when the integral term in equation (8) is 
neglected. We will consider this approximation in the next 
section. For this we will need to have the explicit expression 
for the matrix 7. Let us derive it now. 

The derivation can be based on matrix techniques [4, 6] of 
the solution of the system of ordinary differential equations (5). 
This is a straightforward way. However, it leads to 
quite cumbersome final expressions, if we consider arbitrary 
matrices G,,,; with all non-zero elements. Thus, we use another 
approach, based on the relationship between the components 
of the electric field 


E = Eë + Ena, (9) 


where E; and E are vectors, which are parallel and 
perpendicular to the meridional plane defined by the axis 
OZ and light propagation direction, and the components 
L, h, h, I4 of the Stokes vector [2] I: 


l = E, E} + E2E%, (10) 

h = E, Ej — E2E3, (11) 

h = E1 Eš + E* E, (12) 

l, = i(E1 E3 — EXE), (13) 

where we have omitted a common multiplier. Quadratic 


forms (10)-(13) should be averaged, taking into account 
rapid oscillations of electric field (9) for the finite (and 
large compared to the period of oscillations) time of optical 
measurements. 

Clearly, the matrix Î can be obtained directly from the 
solution [3] of the equation for the coherent field E: 


> 


dE = 
— =-ME, 


JE (14) 


where 


M =27r Nk? F (0), (15) 


N is the number of particles in a unit volume of a random 
medium, k = 27/A,A is the wavelength and F(0) is 


the averaged amplitude scattering matrix in the forward 
direction [2,4]. The solution of equation (14) has the following 
form [4]: 

E=CEo, (16) 
where the matrix C is as given in appendix A. Let us find the 
relationship between the matrix T in equation (6) and matrix 
Cin equation (16). 

For this we will use the density matrix [2, 3] 


p=E@E’, (17) 


where @ means the direct product. It follows from 


equations (16), (17) that 
(18) 
or foe 

ô = ChoC*, (19) 
where ĝo = Éo ® Et . On the other hand, the Stokes vector 


components (10)-(13) can be obtained from the following 
formula: 


I; = Tr(®jĝ), (20) 
where Tr means the trace operation and 

a Ql) 

a=(j 9). (22) 

a=({ 4). 23) 

s= oy (24) 


Aside from the permutation of indices, these matrices are 
usually referred to as Pauli spin matrices. These sigma 
matrices have the following properties: 


Oi = 0], 


A+ A 
0; = 0i, 


(25) 
Tr(G;G;) = 2ôij, 


where i, j, k take the values 2, 3, 4 in cyclic order and ô;; is 
the Kronecker symbol. 


We obtain from equations (19), (20) for components of 
the Stokes vector of a coherent part of the light field: 


(Ie); = Tr(6jCPoC*), (26) 
or 
4 A A 
(i); = } Tr p G ôa (27) 
k=1 
where we have used the expansion 
4 
ĝo = DD (64 Jk). (28) 


k=1 


The value of J% is scalar and, therefore, can be taken out of the 
trace operation. Then it follows from equation (27) that 


A A 


L=TJ, (29) 
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where elements of the matrix T are given by 


Tjk = 4 Tr(ê;Ĉô;Ĉ*). (30) 
Equations (29), (30) can be used for studies of coherent beam 
propagation in a disperse medium. For this one needs to know 
the matrix C, which is presented in appendix A. 

For instance, it follows for media with optically active 
spherical particles after simple algebraic calculations that 


cosh bl 0 0 — sinh bl 
z 0 cosal sinal 0 
a 0 — sinal cosal 0 exp(—el), 
— sinh bl 0 0 cosh bl 


where e = 42 Nk~* Re(F(0)), a = 40 Nk~? Re(F)2(0)), b = 
4r Nk~? Im(F)2(0)). Here Fıı(0) = Fo(0) and F)2(0) = 
—F (0) are elements of the amplitude scattering matrix for 
optically active spherical scatterers [2, 7]. 


3. Diffused light 


3.1. General equations 


Let us now consider the diffused or incoherent part of the light 
field in the framework of the single-scattering approximation. 
Then the integral term in equation (8), which accounts for 
multiple light scattering, can be dropped and we arrive at the 
system of four ordinary nonuniform differential equations: 


dI(x) 
dx 


= —GenI (x) + Wa), (31) 
where . ee 

W(x) = sca T (x) J (32) 
and x = z/ cos V. The elements of the matrix T are given by 
equation (30). The system (31) can be solved analytically by 


standard techniques [4, 6] (see appendix B). The answer is 
I(x) = HYG), (33) 


where the matrix Ê is defined by the following equation: 


A A 


&@=H'6..,H. (34) 
Here 
A, 0 0 O 
s | O Ar. O O 
“10 0 Az 0 65) 
0 0 0 M 
is the diagonalized extinction matrix. Correspondingly, 


Ay, Az, A3, and A; are eigenvalues of the extinction matrix 
Gexr- By definition [6], the columns of the 4 x 4 matrix Ê are 
given by the eigenvectors of the extinction matrix Ĝexr. 

The elements of the vector Y in equation (33) are given 
by (see appendix B) 


VG) = Nis Osca sp Fp) Jk, (36) 


where Ñ = Ĥ-! and the summation on indices s, p and k 


(from 1 to 4) is assumed. It follows for the elements of the 
matrix F' (see appendix B) that 


4 
Fi œ) = 4) WV), (37) 
r=1 
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where 
Wh, = Tr[ô, DG Di], Y2, =Tr[6,Di6.D3], (38) 


Wi, = Tr[ô, ôx DF], Wi, = Tr16)D26,D3] (39) 


and 
Vie) = exp[—A,ux] — exp[—A;x + (A; — Arwa] (40) 
A; = A,u 
where 4 
a cos (41) 
cos Uy 


Matrices D, and D are given in appendix A. It is assumed 
that w = 0 for the radiation propagated downwards (3 < 5) 
and œ = zo sec (Zo is the geometrical thickness of the layer) 
for the radiation propagated upwards (3 > =). Note that the 
value of u is negative for the radiation propagated upwards and 
positive otherwise. 

We have from equation (40) for the reflected light (0 > 4) 
at the upper boundary (x = 0, œ = zo sec #) of a layer that ` 


1 — exp [(A; — A;u)zo sec 0] 
Ai = A,u ` 


It follows for the transmitted light at the bottom of a slab 
(x = £, œ = 0) that 


vi (0) = (42) 


exp|=A, ul] — exp[—A;£] 


VO = A; — A;u 


(43) 
where £ = zg sec ?. 

Summarizing, we underline that for the calculation of the 
Stokes vector of the diffused singly scattered light in a given 
plane-parallel slab we need to find first of all the eigenvalues 
and eigenvectors of the extinction matrix in this medium. 
This can be easily done, using a variety of analytical and 
numerical approaches [6]. The second step is the calculation 
of traces (38), (39). Finally, the Stokes vector is calculated 
using equations (33), (36), (37). 


3.2. Asymmetric isotropic media with spherical particles 


Let us demonstrate the procedure for the special case of 
media with optically active spherical particles. The relation of 
extinction and scattering matrices of media with such particles 
to their refractive indices and sizes is well known [7]. In 
particular, the extinction matrix in this case simplifies to the 
following general form [2]: 


E 0 0 =b 
r 0 € a 0 
| a a (44) 
-b 0 0 e 
with Aj = € — b, Ax = e + ia, Az =e — ia, A4 = € + b and 
1 0 0 -1 
^ 0 -i i 0 
H=|o 11 0 
1 0 0 1 


Let us find traces (38), (39) now. First of all we note that 2 x 2 
matrices D; and D> have the following simplified forms for 
the case under consideration (see appendix A): 
x x a, +Ô: a a ô — ô 
D= p=, y= pj =A 


(45) 
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The evaluation of traces (see equations (38), (39)) gives 


Yi = Yi = Yh = Yh = Yh = Yh = Ya 
Va = Yu = -Yh = Yass] (46) 
and 
W533 = Yh = Wy = -Yh =i. (47) 


Other values of Y”, are equal to zero in the case under study. 
Then it follows from equation (37) that 


Vi + Vi 0o o Vi-V; 
sg © vi+vi -ivi-vi 0 
F=3| 0 i@i-vd  vi+Vi 0 
vi — vi 0 0 vi+vi 
(48) 


Thus, the problem of finding the Stokes vector of the diffused 
light in an asymmetric turbid medium is reduced to the simple 
matrix multiplications (see equations (33), (36)). 


3.3. Symmetric isotropic media 


The assumption of asymmetric medium gives F = VE, where 
E is the unity matrix and (see equation (40)) 
exp[— Aux] — exp[—Ax + (1 — u) Aq] 


LO) SS aM 


ACU —u) 9) 


where A is the extinction coefficient. We have accounted for 
the fact that ê = AE for the symmetric media. Then we have 


I = Voscad, (50) 


which is a familiar result [8]. 
Specifically, we have for the reflected light at the upper 
boundary (x = 0, œ = £): 


(1 — u)A l 


yrf = (51) 


Also, we obtain for the transmitted diffused light at the bottom 
of a slab (L = £, œ = 0): 


exp[ spans] — exp [- asa] 


A= uA 


V” = (52) 
Indices ref and tr indicate reflected or transmitted light 
correspondingly. The angles % and Ŷọ are counted starting 
from the positive direction of the axis O Z (downward vertical). 
Thus, the value of uo = cos vp is always positive for the 
radiation entering a layer from the top. The same applies to the 
value of cos Ŷ in the case of transmitted radiation. However, 
the value of cos Ŷ is always negative for the reflected light. 
Thus, introducing u = |cos ?|, we obtain 


1 — exp[—Azo(4 + +)] 


yrer = Ho Uo, (53) 
(Ho + u)A j 
Az Az 
i exp[— a] = exp[— =] 
yee 2 eg, (54) 
(0 — WA i 


Combining equations (50), (53), we obtain for the reflected 
light 
[= exp[—Azo(+ + +)] = 


I = (40 + WA — Hoscad . (55) 


A similar equation is easily obtained for the transmitted light. 
Let us introduce the phase matrix 


(56) 


where ø, is the total scattering coefficient, and the reflection 
vector function 


Isa 


WF’ 


where uor F is the incident light flux density. Then it follows 
from equations (55), (57) that 


(57) 


F woP 


> 1 1 + 
R = ——_}1—- —|— + — Jo, 58 
Tel | ata lh i on 


where Jo =J /aF, t = Azo is the optical thickness and 
wo = o; /A is the single scattering albedo. Equation (58) was 
reported by, for example, Hansen and Travis [8]. This confirms 
our derivations. 


4. Conclusions 


The problem of polarized coherent and incoherent light 
propagation in symmetric and asymmetric media was 
considered. The results obtained are based on the solution 
of the vector radiative transfer equation in the framework of 
the single-scattering approximation. Formulae derived (see 
equations (29), (30), and (33), (36), (37)) allow for a simple 
calculation of the Stokes vector of coherent and singly scattered 
diffused light in symmetric and asymmetric isotropic disperse 
media with arbitrary differential scattering and extinction 
matrices. The power of the general method presented was 
demonstrated for the specific case of spherical optically active 
particles. We also re-derived equation (58), which is often 
used for the case of symmetric isotropic media [3, 8]. 

Results reported can be used for the interpretation of 
light scattering experiments, dealing with optically active light 
scattering media [9-12]. 
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Appendix A. The matrix C 


The matrix C is expressed through elements of the matrix M 
(see equation (15)) by the following equation [4]: 


C=D e+ Dye, (A.1) 
where 
a 1 1 —Aj2 
Dee o , 
L= T= Anäi a —Aj2A2 
(A.2) 


i 1 Ea =) 
2 = ——_—__ ; 
1— A124721 —Ar 1 
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and 


(A.3) 


yı = 4(Mıı + Mz + Y4M2 Mz; + (My; — My”), (A.4) 
y2 = $(Mi + Mz — V4My2M2) + (Mi, — M22)?°). (A.5) 


The value of / is the length of the beam path. It follows in the 
case of a plane-parallel layer that l = z sec Yo, where z is the 
geometrical depth and Wp is the incidence angle. 

The matrix T in equation (30) can be written in the 
following form, using equation (A. 1): 


4 
T= 5) Viv", (A.6) 
r=1 


where 
Wi, = Tr[ĉ; Dig Dt], We, = Tr[ô; Dô D3], (A.7) 


Wi, = Tr[ô; Dog Dt], Wi, = Tr[ô; D26 D}] (A.8) 


and 
v, = exp(—A,z sec Yo). (A.9) 
Here we have used the fact [4] that 
A= yi +yř, A =y +y% 
: : (A.10) 
A3=y, +”, A= y +y. 


It follows for optically active spherical particles [7] that 
Miı = Mx, Miı2 = — Mhn; and, therefore, 


yı = Mıı +iMi, yı = My, — iM, (A.11) 
Aj = Ar =l, (A.12) 

ž i 6, +ô a à a, — 6. 
B= i=, by = B= =. I) 


Thus, we have for the electric vector E of acoherent part 
of wave propagating in random media with optically active 
spherical particles (see equation (16)) that 


E = 4{(6, + 64) exp[- (Mı +iM12)1] 
iMi)I]}Eo. 


Here Ep is the electric vector of an incident wave. 


+ (ôi 


64) exp[—(M, (A.14) 


Appendix B. The solution of the system of ordinary 
differential equations 


Let us solve equation (31). For this we multiply equation (31) 
by the matrix H7!: 


> dt 18. Fa aia 
H oa -H` el + HW (B.1) 
be 
and introduce a new unknown vector: 
Y-AT (B.2) 
Then it follows that 
dř pe 
— = —êY +f, (B.3) 
dx 


where (see equation (35)) 
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ê = A'6..,H (B.4) 
and B = ÑW, Ñ = A™!. Recall that W = Gyca7'J (see 
equation (32)). Then it follows that B =Ñ Ssa ÎJ i 

Thus, we should solve four independent ordinary 
differential equations (B.3), which can be easily done by 
standard techniques. The answer is 


Yi(x) = / exp(—Aj (x — s))B;(s) ds, (B.5) 
where (see appendix A) 
Bi (s) = Nis Osca,sq Typ (s)Jp, (B.6) 
4 
Typ (8) = > Wiper ($), (B.7) 
r=1 
Uy (s) = exp(—Amus), (B.8) 


i=1,2,3,4,u= ae a = 0 for the radiation propagated 
downwards (30 < =) anda = zo sec Ŷ (zo is the geometrical 
thickness of the layer) for the radiation propagated upwards 
(0 > 5). Here and below, the summation on repeating indices 
is assumed. We used equation (30) and accounted also for 


boundary conditions: 
Y;(0) = 0 atu > 0 
Y;(€) =0 atu <0, 


where £ = zg sec ?. 
Substituting equation (B.6) into (B.5), we obtain 


Y; (x) = NisOsca,sq F} dp: (B.9) 
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where i 
i 1 r yi 
Fi œ) =} WV), (B.10) 
r=1 
and 
i exp[—A,ux] — exp[—A;x + (A; — A,u)a] 
vV; (x) = -s- Laa. 


Ai ms A,;u 
(B.11) 
Finally, the Stokes vector is obtained from equation (B.2): 


I(x) = HY(x) (B.12) 
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